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Abstract 

Here we use perturbation techniques based on the averaging method to investigate 
Rabi oscillations in cw and pulse-driven two-level systems (TLS’s). By going beyond the 
rotating-wave approximation, especihcally to second-order in perturbation, we obtain the 
Bloch-Siegert shift of the TLS resonant frequency, in which the resonant frequency increases 
with the driving field amplitude. This frequency shift implies that short resonant 7r-pulses 
in which the Rabi frequency is approximately 40% or higher of the transition frequency 
do not achieve complete inversion in TLS’s. Hence, guided by analytical results based on 
the averaging technique, we propose two methods for obtaining population inversions in 
the TLS driven by short 7r-pulses: one with chirping and the other with pulse shaping 
and near resonance blue-shifted detuning. Both methods minimize dephasing due to the 
Bloch-Siegert shift, reduce the dependance of the excitation of the TLS on the pulse phase, 
and are very effective in achieving complete population inversions. 
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I. INTRODUCTION 


Here we investigate the coherent, cw or pulsed, excitation of two-level systems 
(TLS’s) in quantum mechanics, under the semiclassical approximation. Our main 
purpose is to develop approximate analytical methods with better precision than the 
rotating-wave approximation (RWA) that will help us in designing more accurate and 
faster coherent control schemes than what is currently available in the literature. 

We use the perturbative method known as the averaging method M to go one 
step beyond RWA. Usually this method is applied to periodically-driven dynamical 
systems with two very different time scales (one fast and the other slow), although, 
as we will see, even when the time scales are not very far apart the method is still 
very useful. Furthermore, this method has the added advantage that it can provide 
efficiently both transient and stationary approximate solutions, unlike the usual time- 
dependent perturbative methods of quantum mechanics. The proposed approach is 
considerably simpler than to implement than Floquet theory ra, which technically 
is not applicable to TLS’s excited by pulses. By implementing our method, we 
obtained the first correction to the resonant frequency of the TLS as the pump 
amplitude is increased, an effect known as the Bloch-Siegert shift [5j. 

Quantum control of population transfer in TLS’s is a topic of intense current 
theoretical and experimental investigations [HHn], mostly with the purpose of imple¬ 
menting quantum gates in various physical systems for quantum computing. Hence, 
it is desirable that laser pulses used in the realization of quantum gates be as short 
in time as possible so that they conform with basic theoretical models of quantum 
computing B3H- Furthermore, such short pulses reduce the adverse effect of decoher¬ 
ence, which then reduce the amount of errors involved in computation with qubits. 
On the other hand, the fields become stronger, the Fourier spectrum widens, and the 
control becomes more difficult, with the breaking up of the RWA and, consequently, 
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of the area theorem P3HI3], which depends on the KWA. In such situations, the ef¬ 
fectiveness of short and high-amplitude 7r-pulses depends on the phase of the pulse as 
well. Hence, in order to diminish such difficulties, we designed 7r-pulses, either with 
chirp in the carrier frequency or with pulse shaping with near resonance blue-shifted 
carrier frequency that minimizes the dephasing effects due to Bloch-Siegert shift. 
Analytical and numerical results of this approach of pulse chirping or shaping are 
presented with 7r-pulse-driven TLS’s that yield considerably deeper inter level popu¬ 
lation transfers and have almost no dependance on the phase of the pulses. This adds 
to recent investigations on quantum control of TLS’s that try to achieve maximum 
population transfer either with chirping, which mostly use linearly-chirped pulses 
[HHI6] or with pulse shaping mt 

This paper is organized as follows. In Sec.[TT| we review the results of RWA for the 
Schrodinger equation for the ac-driven TLS and present a naive perturbation method 


for tentatively going beyond RWA. In Sec. Ill we obtain the 2nd-order averaging 


method correction to the RWA time-evolution. In Sec. IV we obtain the density- 
matrix equations of motion in the 2nd-order averaging approximation. In Sec. |VJ 
we apply our theoretical results in the design of effective 7r-pulses. In Sec. [VTJ we 
discuss possible applications of our main results and draw our conclusions. 


II. RABI OSCILLATIONS IN A TWO-LEVEL SYSTEM 

The problem of coherent control of TLS’s has a long and intense history, dating 
back to the 30’s. It started with the study of the dynamics of magnetization of nuclei 
with spin 1/2 interacting with a magnetic field B(t) = Bd c + B ac cos(u}t ), in which 
there is an angle between the dc and ac components of the magnetic field. This was 
originally investigated and solved in the RWA, in a semiclassical approach, by I. Rabi 
in Ref. [IE]. He found that, in resonance, there is a precession of the magnetization 
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around the ac field, which became known as Rabi cycle, at a frequency g/dBB ac /h, 
known as Rabi frequency, in which g is the gyromagnetic constant and //g is the 
Bohr magneton. This system has been found to have a time evolution equivalent to 
that of a TLS (atom, molecule, quantum dot, exciton, etc) interacting with a single 
mode of the electromagnetic field inside an optical cavity. 

This generic Hamiltonian for this dynamical system can be written as 

n(t) = £i|l)<l| + E 2 12)(2| + 2M2 0 cos (tat) [|1)(2| + |2)(1|]. (1) 

The wave function in the two-level approximation is given by 

|^(t)) = C 1 (t)|l)+C 2 (t)|2). 

Hence, the corresponding Schrodinger equation is 

ihC\ = EiC\ + 2 hfl 0 cos(ut)C 2 , (2a) 

ihC '2 = E 2 C 2 + 2Hn 0 cos(cot)Ci. (2b) 

With the transformations Ci(t) = e~^ t Ci(t) and C 2 (t) = e~~f?' t c 2 (t), we obtain 
the interaction picture dynamics, which is given by 

ic\ = 2fl 0 e~ lUot cos(cot)c 2 = Q 0 e iSt (l + e~ 2lU}t ) c 2 , (3a) 

ic 2 = 2fi 0 e lWot cos(ut)ci = Q 0 e-‘ M (l + e 2lut ) ci, (3b) 

where u>q = E ' 2 ~ El is the transition frequency and 5 = u> — u 0 is the detuning. 

In general this system would be quasiperiodically driven, and thus one could not 

readily apply the averaging method or Floquet theory. One fixes this problem with 
the following transformation 

i8t 


a(t) = e 2 bi(t), 

(4a) 

iSt 


c 2 (t) =e 2 b 2 (t), 

(4b) 
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we obtain the equations of motion 


ib\ — -bi + h2o (l + e 2ia,t ) b- 2 , (5a) 

ib-2 = — -&2 + fio (l + e 2 * wt ) b \. (5b) 


which are a parametrically-driven dynamical system with periodic coefficients. 
Hence, we can now perform the usual RWA and obtain 


ibi — -b\ + ft 0 b 2l 

£ 

ib 2 = — —&2 + 


( 6 ) 

(7) 


For this approximation to be valid, we assume S and = C(e), with 0 < e << 1. 
With a few algebraic steps one finds 


k = 


5 2 ' 

^5 + ) bi, 


where i = 1,2. For the initial values Ci(0) = &i(0) = 1 and C 2 (0) = b 2 (0) 
are equivalent to &i(0) = 1 and &i(0) = —%5 /2 or to & 2 (0) = 0 and & 2 (0) 
have 


0, which 
—A2 0 we 


bi{t) = cos(flf) — — sin (fit), 
b 2 (t) = --^sin(flf), 


where 



Hence the occupation probabilities in each level are 

I Cl (t) 1 2 = cos 2 (fit) + ^2 sin 2 (fit), 

\Ci{t)\ 2 = sin 2 (fit). 


(8a) 

(8b) 


(9a) 

(9b) 
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In Fig. [lj we compare the ground state occupation probability (population) evo¬ 
lution from the numerical integration of Eq.(|2]) with the corresponding population 
evolution given by RWA in Eqs. The numerical integrator used was Python’s 
Scipy odeint with a time step of 2 -10 . One sees that there is considerable dephasing 
between the RWA and numerical integration results. In order to reduce this type 
of problem, one has to go beyond the RWA. One could attempt to achieve that by 
simply replacing the RWA solution, given in Eqs.(|8]), and using the transformation 
Q backwards, in the right hand of Eq. (J 3 J) . We obtain 

Ci(t) = 1 — 2 ifl 0 [ e~ luloS cos(u}s)e~ l5s ^ 2 b 2 (s)ds 


= 1 

= 1 


2^1 

~TT 


:> -i{ui 0 +5/2)s 


cos(o;s) sin(Os)cis 


'0 


02 ft 

1 0 I p -i(uio+8/2)s 

n 


{sin[(a; + 0)s] + sin[(f2 — o;)s]} ds. 


(10a) 

(10b) 

(10c) 


/ 0 


We need to find the integrals 

ft 


e " v ~ u 'sin[(f] + ui)s]ds — 


'0 


3 -i( u 0 +5/2)s , 
D -i(Q.+2u)-8/2)t _ 

20 + 4c o — 5 


e -i(u 0 +5/2)s gin p _ u j s } ds = 
e ~i(Q-5/2)t _ ^ 

20-5 ' 


g«(0+5/2)t _ ^ 

20 + 5 


_ 2u)+8/2)t _|_ 

20 — 4a; + 5 


(ii) 


( 12 ) 


This simple procedure fails to achieve more accurate results than the RWA, though, 
as one can see in Fig. [2j since it does not decrease the dephasing problem of the 
RWA. Another approach has to be devised in order to give a better approximation 
than the RWA to the Rabi oscillations without too much burden. Since the RWA 
is equivalent to the averaging method in the lst-order approximation, a natural 
candidate to achieve better results than the RWA is the averaging method to second 
order. That is what is accomplished in the next section. 
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III. RABI OSCILLATIONS USING THE AVERAGING METHOD TO SEC¬ 
OND ORDER 


Although one could apply Floquet theory to solve the system of differential equa¬ 
tions (J5]) exactly, such as was performed in Refs. OT, here, instead, we will use the 
simpler approach of going just one step further than the RWA in solving this system 
approximately using the averaging method to second order mm- 
Using the notation of Ref. [13], we have 


0 e~ 2lut 


Vi 

e 2iUJt 0 


2/2 


= -A^o 

with the weakly nonlinear transformation b{t) = y(t ) + eW(y, t ), with b(t),y(t ) G C 2 , 
where W(y,t) obeys the differential equation 

W t = f. (13) 


Hence, upon integration of Eq. (13), we obtain 




Ho 

2c j 


0 e 
e 2l0Jt 0 


e~ 2lut y 2 

—e 2iuit yi 


—2iu)t 



2/i 


2/2 


(14) 


in which the constants of integration were set to zero. The Jacobian matrix of / is 

0 —2zcj£ 


eDf = -Ad 0 


0 e 

0 2lU)t Q 


With another weakly nonlinear transformation y{t ) = z(t ) + e 2 v(z,t), we obtain the 
second-order averaged system, corresponding to the original nonautonomous ordi¬ 
nary differential equation (ODE) system of Eq. (J5]) , 

i = ef 0 (z) + e 2 D/VV, 
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where the overline denotes the time averaging. We then find 


e 2 DfW 



-1 0 


Zl 

2 oj 

0 1 


z 2 


Consequently, we find 


Zl 


1 ^0 


Zl 

zQq 

1 

1 

o 


Zl 

Z 2 


O 0 -§ 


z 2 

2 uj 

0 1 


Z 2 



1 


Zl 






O 0 


z 2 





(15) 


where 6 = 6 — One sees that the Eqs.(15) have the same form as the Eqs. (Td), 
with the correction for the detuning. This correction reduces the dephasing of the 
usual RWA for the Rabi oscillations. This shift in resonant frequency relative to the 
RWA resonant frequency is known as the Bloch-Siegert shift [3, 5J. That is not all 


though, the initial values of the ODE system (15) are different from the initial values 
of Eqs. (JtJ) . Hence, that implies that the solution of Eq. (15) is given by 


z\{t) = A\ cos(flf) + B\ sin(fk), (16a) 

z 2 (t ) = A 2 cos(fit) + B 2 sin(fk), (16b) 


where the coefficients can be found from the initial values 


d(o) 


M°) 

Mo) 

M°) 


A 1 = 


l + l ? 5 


Ao — 


2 u 
flo / (2cu 


i +1 — 

' 2u; 


2 ’ 


2 i 


ORi = — 
2 

ClBo = -l 


(521 (0) + 2O()22(0) 

2fi 0 z I (0) — fe 2 (0) 


(17) 


(18) 

(19) 

( 20 ) 































IV. THE DENSITY MATRIX EQUATIONS OF MOTION 


In order to investigate one qubit dynamics with depopulation and decoher¬ 
ence effects, one has to use the density-matrix equations of motion, instead of 
the Schrodinger equation. We start with the following equations 


A = —71 (A — A 0 ) + 41m crF(t), ( 21 a) 

& = — 72 a + iujoct — iAF(t), (21b) 


where F(t) = 2Q 0 cos (cut+ <fio), A is the population difference W(t)— N 2 (t), A 0 is the 
equilibrium population difference, and a(t) is the off-diagonal element p 2 i{t) of the 
density matrix in the basis of the two states |1) and |2). The parameters 71 = 1/Ti, 
the depopulation rate, and 72 = 1/T 2 , the decoherence rate, are dissipation rates 
obtained from the relaxation-time approximation. We now proceed to implement 


the averaging method with the transformation a(t) = e lujt a(t) to set the Eqs. (21) in 
the slowly-varying frame (also known as the interaction representation). Hence, we 
obtain 


A = - 7l (A - A 0 ) - 2A2 0 (ae lut - d*e~ lut ) (e l{ut+(po) + e~ liu]t+<Po) ) , ( 22 a) 

b = - 7 2 d - x5a - A 1 0 A [e^° + e -^ t+(po} ] , ( 22 b) 

= — 7 2 d* + i5a* + itlo A [e~ lV0 + e l(2ujt+ipo) ] , ( 22 c) 

where we included the equation for a*, because it will be used lates in the averaging 
method to second order. Doing lst-order averaging, which is equivalent to the RWA, 
we obtain 


A = - 7l (A - A 0 ) - 2A2 0 [ae-w - a*e lipo ] , 
a = — 72 U — i5a — zQoAe*^ 0 , 
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(23a) 

(23b) 



(24) 


whose fixed-point solutions are given by 

A2 0 Ae vo zh 2 0 Ao (72 — i8)e l ' po 
72 + i5 7 1 + <5 2 + 4fig 72 / 7l ’ 

Aq 


a = 


A = 


1 + 


4^q72 


(25) 


7i(72 +^ 2 ) 

Next, we proceed to the calculations for the 2nd-order averaging corrections in the 
dynamics of the Bloch equations. 


A. Averaging to 2nd-order 


Again using the notation of Ref. [T9] we obtain 

2 (cre4 2 ‘ Jf+¥3 °) — cr*e -l ( 2wt+¥ ’°)) 
ef = —zO 0 Ag-^wt+^o) 

_Ae*( 2wt +v , o) 

with the weakly nonlinear transformation b(t) = y(t ) +eW(?/, i), with b(t),y(t ) G C 3 , 
where W(y,t) obeys the differential equation 


W f = /. 


(26) 


Hence, upon integration of Eq. (|26), we obtain 

2 (<7e^ 2a ' t+¥ ’°) + 5-* e -*( 2wt +¥ , °)^ 
_A e -!(2wi+ifo) 
_Ae4 2wt +v 3 o) 




(27) 


in which the constants of integration were set to zero. The Jacobian matrix of / is 

q 2e*( 2wt + < / , o) _2 e ~ l ( 2wf +¥ , o) 

eDf = -iQ 0 e -i{2ut+v o) o 0 

_ e t(2uit+cp 0 ) Q Q 
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0 


e 2 DfW 


id 2 

OJ 


a 


-<r 


Hence, in the second-order averaging we obtain 


A = - 7 i(A - A 0 ) - 2iQ 0 [be"^ 0 - a*e vpo ] , 


a = — 7 2 a - 1 



a — tfl 0 Ae l,fo . 


(28a) 

(28b) 


One can readily see that in this approximation the driving phase has no effect on 
the population dynamics, and just adds a constant phase to a. 


V. APPLICATION TO QUANTUM CONTROL 


We now proceed to apply the averaging method to study the action of pulses 
on the dynamics of TLS’s. In particular, we are interested in designing 7r-pulses 


to accomplish the deepest population inversions. By examining Eq. (28), one sees 


that the usual area theorem HM3] does not work. A pulse with central frequency 
at co = u 0 loses resonance as the pulse amplitude increases and, consequently, if 
the peak amplitude is large enough, there is substantial detuning and dephasing 
that population inversion can no longer be accomplished. Furthermore, the phase 
of the pulse may become very relevant in determining the effectiveness of the pulse 
excitation, especially so the shorter the pulse is, whereas in the RWA or in the 


averaging approximation to 2nd-order, it plays no role. Examining Eqs. (28) one 


sees that if the detuning of the pulse is zero during the driving, then one could 
achieve population inversion. Based on this theoretical result, one could devise two 
strategies: one with chirped pulses that take into account the amplitude dependent 
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detuning of the Bloch-Siegert shift, and the other one with pulse shaping instead of 
chirping. 

In the first strategy one could have a chirp on the pulse central frequency given by 
uj(t) = ujq + Qo(t) 2 /cu(t), hence the area theorem could be valid again. This should 
be the case, provided the time variation of co(t) is not too fast. Next, we will use 
only the positive root of this quadratic equation, which is given by 

u(t) = m (29) 

Furthermore, one needs h2 0 (£) to vary slowly in time compared with 2h/ujq and 
flo(t) = 0(e) and uj(t) = 0(e) at most, so that the averaging techniques can be 
applied. One should notice that when TLS’s are driven by pulses, Floquet theory is 
not applicable, whereas averaging techniques may still be applied provided the above 
restrictions are maintained and that the transformation to a rotating frame is given 
by 

a(t) = e'fi uWd8 a(t), 


as was discussed in Ref. [20]. In order to achieve maximum population transfer, a 
Gaussian 7r-pulse should be given by 

u>(s)ds + 

(«-*o) 2 

= 2Q 0 e 2ct o cos(cu 0 t + tp(t)), (30) 


F(t) = 2Q 0 (t) cos 


where the pulse instantaneous frequency cu(s) is given by Eq. (29) and ctq = 


7T 


2 y/2fl c 


and to » (Jo¬ 
in Fig. [4j one sees that this strategy yields very good results, with deep pop¬ 
ulation inversions and little dependance on pulse carrier phase. Whereas, 7r-pulse 
excitations with central frequency lu = u 0 yield results that are very dependent on 
pulse phase. This phase dependance becomes a more relevant issue the shorter the 
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pulse is. The drawback of this method is that it is likely not very easy to design 
the chirped pulses experimentally. In Fig. [5j with the presence of decoherence under 
the relaxation-time approximation the proposed strategy of using pulses with chirp 
yield deeper population invertions than what would be achieved with a resonant 7 r- 
pulse excitation. One obtains a very good fit for the population inversion with the 
second-order averaging result and the numerics with pulse chirping. 

The second strategy, though, involves no chirp, instead it keeps the effective detun¬ 
ing as small as possible by choosing the pulse maximum amplitude so that flo = V / &u. 
Here we need a near resonance condition of 5 = 0(e 2 ), so that O 0 = O(e), since 
u = 0(1). Obviously, for this strategy to be valid, one needs a blue-shifted central 
frequency with respect to the transition frequency ujq = 1 (in our units). The advan¬ 
tage of this method is that one can have an effective 7r-pulse even with a blue-shifted 
detuning, what is not possible when the McCall-Hahn’s area theorem is valid, since 
it is based on the RWA. 

In Fig. [6j we show a sequence of three Gaussian 7r-pulses each one with a central 
frequency of lo = 1.1 and amplitude 2O 0 = 0.2 acting on the TLS. Due to the 
detuning almost only about 40% of the ground state population is excited after the 
three 7r-pulses. Both the RWA and averaging to 2nd-order predict well the behavior 
of the population after the first 7r-pulse, but break down afterwards as can be seen 
in frame (a). In frame (b) we show the instantaneous pulse shapes along with their 
envelopes and in frame (c) we show the accumulated area of the pulse envelopes 
as a function of time to attest that indeed we have 7r-pulses. One sees that if the 
pulse shape is not well chosen the 7r-pulses are not effective, the RWA and averaging 
approximation break down, and one loses quantum control over the qubit. 

In Fig. 0 we again show a sequence of three Gaussian 7r-pulses, each one with a 
central frequency of a; = 1.1 and envelope amplitude with fl 0 = \f8uj = 0.331662 
acting on the TLS. One sees that due to the detuning, after three 7r-pulses the RWA 
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(red line) breaks down yielding less than 80% of population inversion. Whereas the 
averaging to 2nd-order still got approximately 98% of population inversion. With 
this value of pulse envelope amplitude, one achieves nearly 99% of population inver¬ 


sion by numerically integrating the density matrix equations of motion in Eqs. (21). 
This indicates that the averaging method to second-order is a good predictor of the 
behavior of the population transfer between the ground and excited levels when our 
pulse shaping strategy is implemented. In frame (b) we show the instantaneous pulse 
shapes along with their envelopes and in frame (c) we show the accumulated area of 
the pulse envelopes as a function of time to attest that indeed we have 7r-pulses. 

In Fig. [8j we show a sequence of three Gaussian 7r-pulses each one with a central 
frequency of oj — 1.1 and amplitude of twice Go = 0.5 acting on the TLS. Due to the 
detuning almost only about 44% of the ground state population is excited after the 
three 7r-pulses. Both the RWA and averaging to 2nd-order break down as can be seen 
in frame (a). In frame (b) we show the instantaneous pulse shapes (with zero phase) 
along with their envelopes and in frame (c) we show the accumulated area of the 
pulse envelopes as a function of time to attest that indeed we have 7r-pulses. Again, 
one sees that if the pulse shape is not well chosen the 7r-pulses are not effective, the 
RWA and averaging approximation break down, and one loses quantum control over 
the qubit. 

The results presented in Fig. [9] show the amount of population remaining in the 
ground state after a Gaussian 7r-pulse excitation is applied to the TLS. The shapes 
of the pulses in frames (a)-(c) are equal to the shapes of the the pulses in Figs. [b]{8j 
respectively. In frame (a) there is no complete population inversion (under 80%) 
and the outcome of the pulse excitation depends considerably on the phase of pulse. 
Whereas, in frame (b), we notice that our pulse shaping strategy to generate the 
deepest inversions works very well. Here, due to the Bloch-Siegert shift, the pulse 
is effectively in resonance, hence it is a 7r-pulse. We obtain very good fit with the 
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averaging and the numerical results with nearly complete inversion independent of 
pulse phase, what make performing quantum control of qubits is easier. In frame (c), 
since one gets again out of resonance due to the Bloch-Siegert shift, not only there 
is no complete population inversion, but there is also strong influence of the phase 
of the pulse on the outcome of the 7r-pulse excitation. 


VI. CONCLUSION 

Here we showed that perturbation techniques based on the averaging method 
can efficiently describe the unitary (or near unitary) time evolution of TLS’s under 
coherent excitation. With the averaging method to second-order approximation, we 
went beyond the RWA in analyzing the dynamics of such a fundamental process in 
quantum mechanics. The numerical time-evolution of the ground state population 
of the TLS under cw driving had a far better fit with the averaging method to 
2nd-order approximation than with the RWA time evolution. Basically, this better 
approximation was achieved by eliminating the dephasing that occurs in the RWA. By 
going to second-order in the averaging method, we obtained the Bloch-Siegert shift of 
the TLS resonant frequency, in which the resonant frequency increases quadratically 
with the driving field amplitude. This frequency shift is the main cause of the 
dephasing of the ground state population time evolution observed in the RWA. 

Although, the proposed method is not exact as Floquet theory is in the case of cw 
driving of TLS’s, it is far easier to implement than the latter theory. Furthermore, 
it can still be applied even when Floquet theory can no longer be applied, such as 
when the TLS is driven by short pulses. 

We verified that the Bloch-Siegert frequency shift that occurs in cw driving, also 
has important influence on pulse-driven TLS’s. This shift basically implies that the 
well-known area theorem of McCall and Hahn breaks down when the Rabi frequency 
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of pulse-driven TLS’s becomes approximately of the order (> 40%) of the transition 
frequency. By numerically integrating the optical-Bloch equations we note that, 
for intense and short 7r-pulse excitations, consistently more population inversion is 
achieved and with less phase dependance if the pulse central frequency has a chirp 
or the pulse amplitude is chosen so as to minimize the detuning effect of the Bloch- 
Siegert shift. 

The present approaches of designing 7r-pulse shapes guided by the averaging 
method could be used in experimentally generating shorter 7r-pulses that achieve 
considerably more population inversion than with 7r-pulses without chirp or with the 
wrong shape. Our results may prove to be very helpful in speeding up control schemes 
of quantum qubits in quantum computation, since it minimizes the effects of pulse 
phase and decoherence on the effectiveness of the 7r-pulses. If one only achieves an 
imperfect population inversion, such as around 95% of inversion, after a few 7r-pulses 
one could lose the control over the qubits. The control of qubits for the realization 
of quantum gates is of fundamental importance in quantum computation [10], as 
they are usually realized by sequences of laser pulses in physical systems, such as, 
for instance, ion traps !2U (22]. It is important to mention that our methods are not 
restricted to the design of just the Gaussian-shaped 7r-pulses shown here, but they 
may be realized with other types of pulses used in the implementation of quantum 
gates. Furthermore, our results may prove to be useful for increasing the photocur¬ 
rent in quantum dot (QD) devices [23] even when the pulse carrier frequency is not 
in resonance with the interlevel transition frequency of the QD. 

Finally, we would like to point out that the fairly straight forward and powerful 
capabilities of the averaging method, such as those used here, should allow it to have 
a place among the usual time-dependent perturbative methods taught in standard 


16 


text books of quantum mechanics. 
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Rabi oscillations, fl 0 =0.1 



FIG. 1. Comparison of the ground state population time evolution as predicted by the 
RWA in Eq. © and the result given by numerical integration of Eqs. ([3]). In our units 
w 0 = 1. 
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FIG. 2. Comparison of the ground state population time evolution as predicted by the 


naive perturbative method of Eq. (10) and the result given by numerical integration of 
Eqs. ([3]). Note that there is still the same problem of dephasing seen in Fig. (jT]) , showing 
that this perturbative approach does not improve on the RWA. 
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Rabi oscillations, fl n =0.1 
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FIG. 3. Comparison of the ground state population time evolution as predicted by the 


2nd-order averaging approximation of Eqs. (16) with the corrections of Eq. (14) and the 
result given by numerical integration of Eqs. ([3]). One can still see some dephasing, but it 
is far less pronounced than in the RWA results of Fig. [l| 
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FIG. 4. Population in the ground level after a 7r-pulse excitation as a function of pulse 
carrier phase. Before the 7r-pulse excitation N\ = 1. The thick blue line depicts the 
response due to a pulse without chirp, with central frequency cj = 1, while the thin red line 


depicts the response to a chirped pulse with the central frequency w(t) given by Eq. (29). 
The dashed black line is the corresponding 2nd-order averaging result, corresponding to 
an ideal population inversion independent of phase. 
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FIG. 5. Inter level transitions driven by 7r-pulses with dissipation rates in the relaxation 
time approximation. In frame (a) we show the population of the ground level as a function 
of time. The thin blue line depicts the numerical result of the TLS excitation due to a pulse 
without chirp, with central frequency u = 1, the thick red line depicts the numerical results 


of the response to a chirped pulse with the central frequency ui(t) given by Eq. (29), and the 
dashed black line represents the corresponding averaged equations to second order with the 


corrections of Eq. (14) adapted to include chirp. In frame (b) we show the instantaneous 


chirped pulse along with the Gaussian envelope. In frame (c) we plot the central frequency 
and the pulse phase angle as a function of time. 
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FIG. 6. Interlevel transitions driven by three 7r-pulses. In frame (a) we show the population 
of the ground level as a function of time. The excitations are pulses with carrier frequency 
at u = 1.1 and peak amplitudes with Ho = 0.1. The width and phase of the pulses are 
indicated in the figure. In frame (b) we show the instantaneous pulse shapes along with 
the 7r-pulse Gaussian envelopes. In frame (c) we plot the accumulated area of the pulse 
envelopes as a function of time, which shows that we actually have 7r-pulses. 
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FIG. 7. Interlevel transitions driven by three 7r-pulses. In frame (a) we show the popula¬ 
tion of the ground level as a function of time. The thick blue line depicts the numerical 
integration result. One sees that complete population inversions are achieved, despite the 
fast dynamics, in which the RWA result (red line) is no longer valid. On the other hand, 
the results of the averaging method to second-order (dashed line) are still very useful in 
designing effective 7r-pulses. In frame (b) we show the instantaneous pulse shapes with 
central frequency at u> = 1.1 along with the 7r-pulse Gaussian envelopes with peak ampli¬ 
tude with Ho = \ZSlo- One sees that one gets about only one and a half cycles per pulse. 
In frame (c) one sees the accumulated area of the pulse envelopes as a function of time, at 
each pulse the area increases by 7 r. 
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FIG. 8. Interlevel transitions driven by 7r-pulses. In frame (a) we show the population of 
the ground level as a function of time. The excitations are pulses with carrier frequency 
at u = 1.1 and peak amplitudes with Ho = 0 . 5 . The width and phase of the pulses are 
indicated in the figure. The RWA and the averaging results both break down here. In frame 
(b) we show the instantaneous pulse shapes along with the 7r-pulse Gaussian envelopes. In 
frame (c) we plot the accumulated area of the pulse envelopes as a function of time, which 
verifies that we actually have 7r-pulses. 
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FIG. 9. Ground level population as a function of pulse phase after various 7r-pulse excita¬ 
tions with carrier frequency at oj = 1.1. Note that before the excitation IVi = 1 and that 
all pulse envelopes are Gaussian with amplitudes and widths indicated in each frame. In 
frame (a) we have dephasing due to the Bloch-Siegert shift and thus the outcome of the 
7r-pulse is not a complete inversion and is very phase dependant. In frame (b) when the 
pulse amplitude is 2flo = r L\fbio ~ 0.6632 the pulse is effectively at resonance, with nearly 
no dephasing due to the Bloch-Siegert shift. One sees that for this value of amplitude, 
one gets nearly complete inversion and very little phase dependance. This indicates that 


the pulse shaping procedure based on Eqs. (28), i.e. by keeping the detuning as low as 

possible, is very effective. In frame (c) we again have dephasing due to the Bloch-Siegert 
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shift and thus the outcome of the 7r-pulse is very phase dependant. 






























































